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Abstract 

We propose a simple method to compute the velocity difference statistics in forced Burgers 
turbulence in any dimension. Within a reasonnable assumption concerning the nucleation and 
coalescence of shocks, we find in particular that the 'left' tail of the distribution decays as an 
inverse square power, which is compatible with numerical data. Our results are compared to those 
of various recent approaches: instantons, operator product expansion, replicas. 
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1 Introduction 



Burgers equation, which describes the potential flow of a fluid without pressure, provides a wonderful 
laboratory for testing new ideas and techniques in view of the study of fully developped turbulence 
in the Navier Stokes equation. These are two cases of non linear stochastic equations which share 
the same structure of the non linearity. The important difference comes from the nature of the large 
scale structures. In the case of Burgers equation these are shock waves and the corresponding physical 
picture of the flow is rather simple. This simplicity has already allowed for a very detailed study of 
the decaying turbulence jlj [2| . 

The forced case, in which the fluid is stirred randomly and steadily on large length scales, is more 
complicated. However it has been attacked recently by various methods like the operator product 
expansion []|], direct probabilistic methods gj, instanton calculus || and replica method |7|. The 
latter method allowed to get a detailed solution in infinite dimension, and the finite dimensional 
solution seems to be within reach. This would be an important milestone for several reasons. It 
gives an example of a flow with strong intermittency, created by large scale structures. It provides a 
benchmark to test new -or older- ideas on fully developped turbulence. Furthermore this problem is 
also related to interesting problems in condensed matter physics, like the elastic lines in random media 
(e.g. vortices in superconductors), and growth problems ||. In this respect it is interesting to notice 
that Burgers equation (with time playing the role of the running length scale) also appears naturally 
in the renormalization group study of manifolds in random media ||. Therefore, the phenomenology 
of Burgers equation might be directly relevant to experimental studies of pinned Bloch walls |jll| (or 
other elastic manifolds), besides its more direct potential applications to turbulence in one dimensional 
fluid flows [[jl]] or to pattern formation in astrophysics JT2"| . 

In this paper we shall focus on a simple aspect of forced Burgers turbulence: the tails of the velocity 
gradient distribution in the regions where there are no shocks. We evaluate the 'right' tail through a 
rather simple computation, and compare it to the more sophisticated approaches developped recently 
[0, ||, |^|. We then give a conjecture on the 'left' tail which is based on a plausible argument, requiring 
the system to reach a stationnary state. We shall first discuss the one dimensional case, then turn to 
higher dimensions, and compare our results with the previously available ones. 



2 The slope dynamics in one dimension 

In one dimension we consider a velocity field v(x, t) governed by Burgers equation: 

dv dv d 2 v „. . . . 

where f(x,t) is a random force, which is supposed to have Gaussian distribution, with zero mean and 
a second moment: 

(/(*, t)f(x', 0> = eS(t - t')R {^^~\ (2) 

where R is any smooth function decaying to zero fast enough at large arguments (e.g. an exponential 
as in 0]), A is the length scale of the stirring force and e is the injected energy density. To keep 
consistency with the notations of our previous work |7]] , we fix the normalisations by requiring that at 
short distance R{y) = 1 — |y + 0{y 2 ). 

Our approach is based on the following observation: even in the forced case, the flow is organised 
in some smooth regions separated by shocks. Inside the smooth regions, the flow locally depends 
linearly on the position: 

v(x,t) ~ \(t)(x - x Q (t)) . (3) 
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A linear expansion of the original equation allows to derive the evolution of xo(t) (which is irrelevant 
for our present discussion) and of the local slope X(t). This slope follows a Langevin equation: 

^ = -X\t) + r ] (t) (4) 
with a noise term rj, due to the random stirring force, which has a Gaussian distribution, with: 

WW)) = |^(* - • (5) 

This Langevin equation describes the relaxation of a particle at temperature T = in a potential 
V(A) = A 3 /3. The zero temperature case corresponds to the decaying Burgers turbulence, where the 
slope decays in time as X(t) = Ao/(l + Xot), leading to an asymptotic behaviour in which the slope 
is 1/t independently of the initial conditions jjj]. The forced case T / leads to runaway solutions: 
starting from a positive slope Ao, it will be eventually driven by the forcing to an infinitely negative 
slope, A — > — oo. This effect is nothing but the building up of a new shock, which tends to develop as 
soon as the fluctuations due to the forcing drive the system to a negative slope. However the tail of 
the slope distribution at positive A is unaffected by this effect (see below, Eq.(||)): it depends on rare 
noise configurations which, for a Gaussian forcing, lead to a Boltzmann form: 

P(A) ~ Cfe-W = Cexp , (A - +oo) , (6) 

where C is a constant. A similar result was obtained in the context of pinned Charge Density Waves 
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Inside the linear regions, the velocity difference between two points at distance r,u = v(x+r)—v(x), 
equals the slope A times the distance r, so that the slope fluctuations induce the following tail for the 
pdf of the velocity difference: 

m a C ( 2u 3 A 2 \ 

P ( u ) - ~ ex P I — g^r I . (u^+oo) (7) 

which is precisely the result derived in 0]. We shall discuss the relationship of this elementary 
derivation with these other more elaborate ones later on. 

From the evolution of the slope (||) one can also discuss, in a more speculative way, the other 
tail of the pdf of the velocity differences (at large negative u). Since the regime we want to study is 
stationnary, the average number of shocks (or of linear regions) must be conserved. Hence, the slopes 
which 'disappear' at A = — oo creating a new shock and thus two new slopes, must be compensated 
by the spontaneous disappearance of some other slopes through shock coalescence. In other words, 
for the number of 'cells' (locally linear regions) to be conserved, one must introduce a source term 
JqS(X) in the corresponding Fokker-Planck equation, which describes these 'reinjection'/coalescence 
processes. These processes (and hence the correct form of 5(A)) are difficult to describe properly, but 
fortunately the asymptotic results are rather insensitive to the detailed shape of S(X). We normalize 
the integral over A of S(X) to 1, and write the Fokker Planck equation for the Langevin process (Q) in 
the presence of the source term as: 

This equation now possesses a stationnary solution with a non vanishing current Jo, which is nothing 
but the probability per unit time for a new shock to appear. This stationnary solution is given by: 

P(A) = Jo exp (-^\ J" dX'F(X') exp (+^\ (9) 
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where F(X) = f£° dX' S(X') and Jo is fixed by the normalisation of P(X). If we make the reasonable 
assumption that the reinjection term S(X) is a rapidly decaying function, the integral appearing in 
the above form (g) of the slope probability converges for A — > +00, leading back the above result, Eq. 
(^). For A — > —00, F(X) — > 1, and the asymptotic form of P (A) is easily shown to be: 

P(A)~§, (A — > -00) . (10) 
This in turn leads to a velocity difference distribution which decays as: 

P(«)~^, («--oo). (11) 

which happens to be precisely the form that Chekhlov and Yakhot have proposed using very different 
arguments, and which is compatible with numerical data jy|]. 

Let us emphasize that the Fokker-Planck equation without the source term S(X) does not have a 
(normalizable) equilibrium state. This normalizable stationnary state can only appear in the presence 
of a non-zero current for large negative A. The absence of an equilibrium state is also at the heart of 
the operator product expansion approach of Polyakov on which we shall come back below. 



3 The slope dynamics in higher dimension 

The generalization of this approach to the TV-dimensional Burgers equation is interesting. The velocity 
field verifies the equation: 

dv _ ^ ^ 2 



at 



+ (v-V)v = vV z v + f(x, t) 



(12) 



where v is a gradient flow and the forcing term / is a random gradient field, with a Gaussian distri- 
bution of mean zero and second moment given by: 



p(x,t)f k (x',t') = e5(t-t') 



k {x- x'y (x - x 1 ) 



Afe 



NA 2 



G 



(x — x') 



•l\2 



2NA 2 



(13) 



where the N dependance is chosen such as to insure the existence of a well defined large N limit 0. As 
before, the correlation G decreases fast enough at large arguments and behaves as G(y) ~ 1 — | + 0(y 2 ) 
at small arguments. 

Between the shocks (which now have a dimension N — i), the flow is locally radial, with a slope 
matrix M: 

Vj (x,t) =J2 M jk(t)(x k -x° k (t)) . (14) 



The slope matrix evolves in time as: 



dMjk 
dt 



(M 2 ) . k {t)+y 3k {t) , 



where y-f. is a Gaussian random noise with a zero mean and a variance given by: 

(yij(t)yki(t')) = j^i^ij^M + s ik 8ji + Su8jk)S(t - 1') . 



(15) 



(16) 



We shall need to study the statistics of the eigenvalues Aj of the slope matrix M . From second order 
perturbation theory, one finds that these eigenvalues verify a Langevin type equation: 

dt 



-K + 



1 



iVA 2 .f- A; - A, 



+ Vi > 



(17) 
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where the level repulsion term has appeared ||15|| , and the noise is Gaussian with a correlation: 

fa(%(0> = jf&( 1 + 2 %) *(* - • ( 18 ) 

Compared to the one dimensional case, the situation is more complicated because of two effects: level 
repulsion on one hand, and correlation of the noises rji acting on each eigenvalue on the other hand. 
This noise correlation prevents the existence of an equilibrium (currentless) Boltzmann distribution, 
which one could use as above to obtain the right tail of the slope distribution. We shall instead use 
an instanton computation, which is well suited to deal with such correlations (see e.g. [16]). 
We want to compute the pdf P(u) of the longitudinal velocity difference, 



u = ^2 ( v i(x + r) - Vi(x)) -i- , (19) 
i ' ' 

which is given between the shocks by: 



u = rJ2riM ij f j , (20) 

ij 



where r = \r\ and f is a unit vector in the direction of r. In terms of the eigenvalues of M, it takes 
the form: 



N 

u = r 

i=l 



Em 2 ( 21 ) 



(22) 



where a is a random unit vector (obtained from r through a random rotation). 

Let us introduce the Laplace transform of the pdf of the longitudinal velocity difference: 

G(/i,r) = f duP{u)e^ u = f °° ^ f d ai ...da N (exp I fir E a%\j + <E a ) ~ X ) ] > 

J -^ 2mJ \ 3 J I 

In the spirit of recent works in turbulence [17, 18, ^| and in other fields [|^, IS], we shall use a path 
integral representation for the probability distribution of the eigenvalues, and evaluate the large \i tail 
of G(/i, r) by finding the leading instanton trajectory. The path integral is written through standard 
manipulations |2(J. The dynamics is started at a time to ~~ > — oo, and the eigenvalues are measured at 
t = 0. We introduce conjugate fields \i(t) to implement the Langevin equation, and get after average 
over the thermal noise: 

/ioo Ay r ^ r _ 

_ ioo m J II da J J II d[XM*i\ ex P (23) 



with an action: 



+ f dt w& ( 2 E^(*) 2 +E^(*)^(*))+^E4 A i( t = )+ 2 (E a2 - 1 ]( 24 ) 

00 \ 3 j,k / j \ 3 J 

The stationnarity conditions of this action are easily written. The effect of the parameter [i on the 
trajectories Xj(t),Xj(t) is just to impose the boundary condition Xj(t = - ) = —firaj. The action 

S has the following scaling property: under the change fi — > 7 2 /i, t — ► jt, A« — > 7A1, X\ —>■ 7 2 Aj and 
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z — > 7 3 z, all the terms scale as 7 s , except the level repulsion term which is left invariant. Therefore, 
in the large /x limit, it is legitimate to neglect this level repulsion effect, and the following instanton 
solves these simplified equations: 

(22 = ... = fljV = A2 = ... = \n = 

w = t*# ^ t) = -jr^w (25) 

where the parameter <f> is: <fi = y / 3efir/(2NA 2 ). The action of this instanton is: 



- In(G(M, r))*S = -J ^2 (^) 3/2 (26) 
Assuming that this instanton gives the leading contribution to the Laplace transform G at large /j,, 



one deduces from (26) the following tail of the pdf of velocity differences at large u: 



( 2NA 2 u 3 \ 

P(u)~Cexpl — (u^+00). (27) 

In order to check the dominance of the instanton solution (f25|), we have performed two checks: an 
evaluation of the contribution of other instantons, and a numerical simulation of the Langevin equation. 
Without level repulsion, other solutions can be found, with (N —p) a,j and A, identically zero, and the 
p remaining A and A all equal. The action of these instantons is a factor f p = J smaller than (|26|), 
Since f p <l for all p > 1, these other solutions can all be neglected. Taking into account the level 
repulsion can only increase further the action of these other solutions, so they should not contribute. 

A direct numerical test of the result is easily done. We have simulated the Langevin evolution 
equation for the eigenvalues (|l7|) in the case N = 2. We adopted a simple Ito discretization of time 
and computed the histogram oiu/r = X\ cos 2 (ip) + A2 sin 2 (^) where ip is a random angle with uniform 
distribution in [0, 2ir]. The time step has been kept adaptive in order to be able to use a small time step 
whenever the two eigenvalues are close to each other. Our simulation has been done at e/A 2 = 10/9 
and the statistics is over 10 6 points. The result of the histogram is given in Fig.l, where the straight 
line gives the slope predicted from our instanton computation (^7]). The agreement confirms that the 



instanton (25) indeed gives the leading contribution. We have also simulated the Langevin process 
without the level repulsion term, and checked that we find the same tail for lnP(u). Let us just 
mention that the finite u effects turn out to be larger in this latter case, and the analytic prediction 
for the tail is found only in the regime where lnP(u) < —10. This finite u corrections should also be 
important for the analysis of velocity difference tails in numerical simulations of Burgers equation. 

Altogether, we are confident that the right tail of the velocity difference is given by (^7j) also in the 
multidimensional case. It is important to notice that this tail gets thinner and thinner as N grows, 
and finally disappears in the limit N — > 00, which was actually the situation that we considered in our 
previous paper [0] using the replica method. Indeed, the calculation that we performed there shows 
that, for iV = 00, P(u) is strictly zero when u > (where ua = (eA) 1 / 3 is the velocity at scale 



A). This is compatible with Eq. (27) for N = 00. It however clearly emphasizes the limitations of a 
large N approach: although many features of the velocity field are correctly predicted, such a method 
is unable to predict the tails of the velocity distribution. Stated differently, the limits of large velocity 
difference and large dimension do not commute. 

4 Comparison with the direct instanton approach 

The previous results (|6|,^7|) on the tail of the longitudinal velocity difference at large positive u are 
identical to those obtained in H ||] in one dimension and in || in dimensions N > 1. The one 
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dimensional result is thus confirmed by a very simple and direct derivation. In higher dimensions, 
we are unfortunately not able to solve fully for the probability distribution of the ./V eigenvalues 
because of the noise correlations. We have rather found a simple instanton solution, which allows us 
to generalize the one dimensional result to arbitrary dimension. A naive generalisation of Eq. (|8|) to 
higher dimensions also suggest that the left tail of the distribution of velocity difference should decay 
as u~ 2 in any finite dimension. 

Although our instanton computation looks much easier than the one developped by Gurarie and 
Migdal H (we are dealing with the probability distribution of N eigenvalues Aj instead of a N com- 
ponent velocity field Vi(x)), the two methods are actually very similar. In one dimension for instance, 
the instanton for the velocity field found in Q corresponds precisely to a velocity field growing linearly 
with distance, with a slope called a in ||. It is easy to deduce from the equations (19,20) of |6| 
evolution equation for a which is identical to the one which is derived from the stationnarity of 
The same is true in larger dimensions. 

To conclude our comparison with the work of Q, we have found that the positive u tail of P(u) de- 
rived from their nice instanton solution is certainly correct at least in one dimension. Our computation 
also indicates that this instanton cannot directly give the tail at negative u. In fact its structure, with 
a linear velocity field, shows that this instanton takes into account basically the linear regions, and 
not the shocks. As we argued before, these linear regions dominate the right tail of P(u), but the left 
tail is given by a completely different process. In order to understand this left tail quantitatively, one 
needs to control both the formation of shocks, and the statistics of the sizes of the jumps in velocity at 
the shocks (or equivalently the statistics of the lengths of the linear regions) . This has been achieved 
in the decaying problem O, or in the forced case in large dimensions [p], and it is hard to think that 
this information is contained in a simple linear instanton configuration. 




5 Comparison with the operator product expansion 

Let us now discuss the operator product expansion (OPE) approach of Polyakov ||, which is supposed 
to give the full pdf P (u) in the scaling region. Using some conjecture on the structure of the fusion 
rules in the OPE, one finds that, in one dimension, this P(u) is exactly the totally assymmetric Levy 
stable distribution of index |, given by the inverse Laplace transform of the function G(fi, r) found in 
(^6|). The behaviour of the right tail of P(u) is again given by Eq. @, and is thus compatible with 
the other approaches. However one should emphasize that this tail does not constitute a test for the 
fusion rules of the OPE. In the simple case of the 'two point' function P(u), the OPE approach relies 
on the following steps. Restricting here to the case N = 1, one first introduces a generating functional 
Z(si,S2,t) = (exp[siv(xi, t) + S2v(x2, t)]). In the stationnary regime, one should have = 0, which, 
using the tricks of ref.||, can also be written as: 

+ t; t\ z o — = ([*i/(si> *) + s 2 f{x 2 , t)\ exp[ Sl v( Xl ,t) + s 2 v(x 2 ,t)}) + A , (28) 



dsidxi ds2ds2 sidxi S2&X2 

where A is a term coming from the viscosity contribution, which is singular in the v — > limit (the 
"anomaly"). Using the fact that the force / is Gaussian, the first term in the right hand side can be 
reexpressed as: 

^[R(0)(sl + si) + 2R( X1 - x 2 )s lS 2] (29) 



Now, changing variables to: 



Xl = x + - x 2 = x-- [i= — - — fj, = — - — (30) 
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and using translational invariance, one finds, for [/,' = 0, the following final equation for G(/i, r) = 
(exp^d) -«(-§)]): 

where r « A is assumed. The simplest fusion rule leads to a contribution scaling as 

i~aG+^ (32) 
\i or 

where a and (5 are two constants to be determined. It turns out that the right hand tail of the solution 
G(fj,,r) of (31) has a leading behaviour InG ~ ^/2e/(3A 2 )(^r) 3 / 2 which is independent both of a and 



(3 ! Therefore the exp(— cu 3 ) tail comes out exactly as what we derived in (|7|) (with precisely the 
same c), independently of the fusion rule. This is actually expected since we argued that this tail 
comes from the linear regions where one can completely neglect the anomaly. Without anomaly, the 
above equation has no solution which can be interpreted as the Laplace transform of a probability 
distribution (cf. the discussion after Eq. (|i~l])). The fusion rules thus become crucial when one looks 
at the negative u tail. If one insists that the solution G(/z, r) of ([H]) is the Laplace transform of a 
(normalizable) probability law, with an anomaly given by (|32| ) with a = for all fj,, one finds that P 
is a Levy distribution of index |, with a left tail decaying as P(u) ~ |ii|~ 5 / 2 Q. This result disagrees 
with our above calculation, which suggested a \u\~ 2 tail from the existence of a finite outgoing current 
at large (negative) slopes, related to the formation of new shocks. We do not have enough control 
on this 'reinjection process' of the slopes to make any strong statement. However, our argument is 



suggestive enough to try to look at other fusion rules which indeed lead to such a tail [21]. 

6 Conclusion 

We have discussed the tails of the velocity difference distribution in Burgers' forced turbulence. We 
have shown how a direct method gives result which are in agreement with other methods in the right 
tail of the distribution. This method however suggests that the left tail is governed by the dynamics 
of shock formation, and requires some control over the shocks nucleation/coalescence processes. A 
simple hypothesis describing these processes in terms of an effective source term in the corresponding 
Fokker-Planck equation gives a \u\~ 2 tail for the left side of the distribution, at variance with the 
simplest OPE conjecture, but apparently in agreement with numerical simulations. 

It would be very interesting to extend these methods to describe more complicated quantities, such 
as the statistics of the velocity potential h defined as v = — V/i. In particular, the statistics of the 
'barrier heights' separating the different valleys of h (corresponding to the cells of the Burgers flow) 
would be of direct interest to understand in detail the long time dynamics of randomly pinned objects^]. 
In particular, it would be important to see to what extent the large N result, which states that these 
barriers are exponentially distributed, has to be modified for finite N. Another interesting perspective 
is the generalization of these approaches to the case of a forcing correlated on long distances, where 



Kolmogorov's scaling is observed between the shocks [22, p3| 



Acknowledgments. We want to thank A. Chekhlov, T. Dombre, V. Gurarie and V. Yakhot for 
interesting discussions. 

1 Conversely, as originally noticed by Feigel'man |Jl3|, the very large 'slope' tail A — > oo (Eq. (^)) of the Burgers flow 
is related to the high frequency response of these pinned objects 
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Figure Captions 
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1 The logarithm of the probability distribution function of velocity differences, lnP(u), versus 
(u/r) 3 , in dimension N = 2, for e/A 2 = 10/9. The curve is obtained from the simulation of the 
Langevin equation on the slope matrix eigenvalues, and the straight line indicates the prediction 
from the instanton computation (ETj). 
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